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Abstract. The contribution of ordered rotation to the observed tilt and thickness of the Fundamental Plane of 
elliptical galaxies is studied by means of oblate, two-integrals cuspy galaxy models with adjustable flattening, 
variable amount of ordered rotational support, and possible presence of a dark matter halo and of a central super- 
massive black hole. We find that, when restricting the measure of the velocity dispersion to the central galactic 
regions, rotation has a negligible effect, and so cannot be responsible of the observed tilt. However, streaming 
velocity effect can be significant when observing small and rotationally supported galaxies through large (relative) 
aperture (as for example in Fundamental Plane studies at high redshift), and can lead to unrealistically low mass- 
to-light ratios. The effect of a central supermassive black hole on the kinematical fields, and the models position 
in the t>/<T-ellipticity plane are also discussed. 
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1. Introduction 

In the observational three-dimensional space of central 
velocity dispersion co, (circularized) effective radius R e , 
and mean surface brightness within the effective radius 
I c = L/(2irRc) (where L is the total galaxy luminosity), 
early-type galaxies approximately locate on a plane, called 
the Fundamental Plane (hereafter FP; Dressier et al. 1987; 
Djorgovski & Davis 1987), and represented by the best-fit 
relation: 



log R c — a log er + b log I e 



(1) 



The coefficients a, b, and c depend slightly on the con- 
sidered photometric band (e.g., Pahre, de Carvalho, & 
Djorgovski, 1998; Scodeggio et al. 1998): for example, by 
measuring R c in kpc, do m krn/s, and I c in L Q /pc 2 , re- 
ported values in the Gunn r band are a — 1.24 ± 0.07, 
b = -0.82 ± 0.02, c = 0.182 1 (J0rgensen, Franx, & 
Kjasrgaard 1996, hereafter JFK). One of the most striking 
observational properties of the FP is its small and nearly 
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1 This value of c refers to the Coma cluster and to Ho = 50 
km s _1 Mpc -1 . Note that erg is usually corrected to a circular 
aperture with diameter 1.19 h~ x kpc, corresponding to a radial 
range ~ 0.05 R a -R c for h = 0.5, and for typical values of R c 
(e.g., J0rgensen, Franx, & Kjaergaard 1995). 



constant thickness: the distribution of log R c around the 
best-fit FP (at fixed Op an d I e ) has a measured rms (here- 
after cri n t) corresponding to a scatter in R c ~ 15% -j- 20% 
(after correction for measurement errors, see, e.g., Faber 
et al. 1987; JFK). 

In addition, for a stationary stellar system the scalar 
virial theorem can be written as 



Rn 



= Kv a. 



(2) 



where Y„ is the stellar mass-to-light ratio in the photo- 
metric band used for the determination of L and R e , while 
the coefficient Ky takes into account projection effects, 
the galaxy stellar density distribution, the stellar orbital 
distribution (such as velocity dispersion anisotropy and 
rotational support), and the effects related to the pres- 
ence of dark matter. Equations Q and @ imply that 
in real galaxies, no matter how complex their structure 
is, Y*/Ky is a well-defined function of any two of the 
three observables (L, R e ,ao) and this dependence is com- 
monly referred as the "FP tilt". In other words, the FP 
tilt indicates that structural/ dynamical (Ky) and stellar 
population (T») properties in real galaxies are strictly 
connected, possibly as a consequence of their formation 
process. Overall, the interpretation of the FP cannot be 
limited to the study of its tilt only, but requires to take 
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consistently into account also its thinness (as done, for 
example, in the "orthogonal exploration" approach, see 
Renzini & Ciotti 1993, Ciotti, Lanzoni & Renzini 1996, 
Ciotti & Lanzoni 1997). A step forward in the study 
of the structural and dynamical implications of the tilt 
and thinness of the FP was made by Bertin, Ciotti & 
Del Principe (2002), who introduced, although limited 
to spherical models, the more general Monte-Carlo ap- 
proach. The Monte-Carlo approach was successively ap- 
plied by Lanzoni & Ciotti (2003, hereafter LC) to the 
investigation of projection effects on the FP thickness, 
by adopting fully analytical axisymmetric one-component 
Ferrers galaxy models (Ferrers 1877). In particular, LC 
found that while projection effects do contribute to the 
observed FP scatter, nonetheless the FP physical scatter 
(as determined by variations of the physical properties 
from galaxy to galaxy), sums up to 90% of the intrinsic 
FP scatter. In addition, LC found that the contribution 
of ordered streaming motions to the observed velocity dis- 
persion is negligible when small/medium apertures (< R e ) 
are used for the spectroscopic observations. Unfortunately, 
while fully reliable about projection effects on R e (which 
are independent of the specific homeoidal density profile 
adopted) , the central regions of Ferrers models are unre- 
alistically flat, so that doubts can be risen about the use 
of their central velocity dispersion to estimate streaming 
velocity effects. Moreover, what is the additional effect of 
a central super-massive black hole (SMBH) or of a dark 
matter halo? An attempt to answer the questions above 
was made by Ciotti & Bertin (2005, hereafter CB) by us- 
ing fully analytical cuspy galaxy models obtained from a 
homeoidal expansion, qualitatively confirming the results 
of LC. Here, in order to extend the results of LC and CB, 
we numerically solve the Jeans and projected dynamics 
equations for more realistic, oblate galaxy models with a 
central cusp, and variable amount of internal streaming 
velocity, allowing for the presence of a central SMBH and 
a dark matter halo. 

The paper is organized as follows. In Sect. |21 we de- 
rive the relevant intrinsic and projected properties of the 
adopted models, while in Sects. 01 ^ and [21 we present the 
main results of our analysis together with some observa- 
tionally related consequences, such as the model position 
in the v/a-ellipticity plane. In Sect. we summarize the 
results, while in the Appendix we present the fully ana- 
lytical model adopted to test the numerical code. 



2. The models 

2.1. 3D quantities 

At variance with LC, who adopted the centrally flat 
and spatially truncated Ferrers model, we now use cuspy 
oblate galaxy models with homeoidal density distribution, 
belonging to the family of the so-called 7-modcls (Dchnen 



1993; Dehnen & Gerhard 1994, hereafter DG; Tremaine 
et al. 1994; Qian et al. 1995). Their density profile is 



P* = 



47Ti?3, 



(3-7) 

^(l+m) 4 -v (°^< 3 )< 



(3) 



where is the total mass, R c the characteristic scale, 

and in cylindrical coordinates 2 m = \J R 2 + z 2 /q 2 , with 

R = R/R c and z = z/R c ; the parameter < q < 1 mea- 
sures the intrinsic model flattening. Note that these mod- 
els for 7 ~ 0.75 and 7 ~ 1 provide a good representation 
(in the central regions) of the de Vaucouleurs (1948) and 
of the Navarro, Frenk, & White (1997) profiles, respec- 
tively. In two-component models the dark matter halo is 
also described by a density profile as in Eq. © , in general 
with a different total mass, scale-radius, density slope and 
flattening; a SMBH of mass Mbh can be also added at the 
model center. In all cases, the stellar mass-to-light ratio 
T* is assumed to be constant within each model. 

We assume that the density profile in Eq. © is sup- 
ported by a two- integrals distribution function f(E, L z ) 3 
(where E and L z are the energy and the z-component 
of the angular momentum of stars per unit mass, respec- 
tively), and the Jeans equations are 



dp*cr R 

dz 
dp*o 2 R 

dR 



P* o J 

oz 
P*( v l ~ 



4) 



R 



= -P* 



dR' 



(4) 
(5) 



where (f> is the total gravitational potential, v R = v z = 0, 
the off-diagonal elements of the velocity dispersion tensor 
vanish, and o 2 R — o 2 (see, e.g., Binney & Tremaine 1987, 
hereafter BT). 



As in LC, v 2 is splitted into streaming motion v, 



and azimuthal velocity dispersion cr 2 = (v v 
v 2 with the Satoh (1980) /c-decomposition: 



)2 = v 2 - 



= fc 2 (v 



and 
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rl), 



fc 2 )(v 2 -4), 



(6) 



(7) 



with < k < 1. For k = no ordered motions are present, 
and the velocity dispersion tensor is maximally tangen- 
tially anisotropic, while for k — 1 the velocity dispersion 
tensor is isotropic, and the galaxy flattening is due to az- 
imuthal streaming velocity (the so-called "isotropic rota- 
tor" case). 



2 The (R, (p, z) coordinates are related to the natural 
Cartesian coordinates by the relations R = -U 1 x 2 + y 2 , cos ip = 
x/R, sin^ = y/R. From now on, the symbol "~" over a coor- 
dinate will indicate normalization to R c . 

3 Note that the important issue of the models phase-space 
consistency is beyond the tasks of this work (see Qian et al. 
1995). We use symbol v for the velocity in the phase space, 
while v (x) = v is the streaming velocity. In general, a bar over 
a quantity means average over phase-space velocities. 
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As well known, the gravitational potential of the den- 
sity distribution in Eq. (0 cannot be expressed in terms 
of elementary functions, and so we compute it by using 
an expansion in orthogonal function: Eqs. (0J) and © are 
then integrated numerically fSect. EO]! . 

2.2. Projected quantities 

In order to project the model properties on the plane of 
the sky (the projection plane), we define the observer sys- 
tem as a Cartesian coordinate system (xl, yr, zi) with the 
line-of-sight (los) parallel to the zi axis, and with the xi 
axis coincident with the x axis of the natural Cartesian 
system introduced in Sect. 12.11 The angle between z and 
z/isO < 9 < 7r/2: 9 = corresponds to the face-on 
view of the galaxy, while 9 = n/2 to the edge-on view. 
With this choice, the projection plane is (xf yi) and the 
los direction in the natural coordinate system is given by 
n = (0, — sin 6, cos#) 4 . Following LC we indicate with T,(£) 
the projected (mass or light) density distribution, where 



xf 



q(0) 2 



y' 2 /q(d) 

,o2 



q 2 sin 2 I 



(8) 



The quantity £ determines the size of the elliptic isophotes, 
and their ellipticity is given by e = 1 — q(9). For fixed £ 
the major and minor isophotal semi-axes are a = £ and 
b = q(9)£, and the associated circularized radius is de- 
fined as ttR 2 = irab, i.e., Ri = \J q{9)£. In particular, the 
circularized effective radius is 



(9) 



where M p (4) = M/2, with M p (£) being the projected 
mass within £. Thus the identity R c = y/ q(9)R® (where 
R® the effective radius of the model when seen face-on or 
in case of spherical symmetry) is a common property of 
all axisymmetric homeoidal distributions, independently 
of their specific density profile, and so the dependence on 
9 of the circularized effective radius of the present models 
and those in LC is identical. 

The projected velocity fields at (xf, yl) are obtained 
by numerical integration of the projection along n of their 
spatial counterpart. This is done by transforming the cor- 
responding spatial velocity moments from cylindrical to 
Cartesian coordinates (LC, Eqs. [21]-[27]). Note that, in 
presence of a non-zero projected streaming velocity field 
Up, the velocity dispersion accessible to observation is 



' los 



(xf, yf) 



1 



p ((v, n) - Up) dzi = 



V 2 -v 2 



(10) 



where (,) is the standard inner product, and for the def- 
inition of V 2 (the los integration of v 2 , i.e. the squared 

4 The los vector points toward the observer, and so positive 
velocities correspond to a blue-shift; the general framework for 
the projection is given in LC. 



component along n of the streaming velocity field) and of 
a 2 , we again refer to LC. Note also that, independently of 
the los orientation, v n = on the isophotal minor axis yi, 

a 2 v every- 



and thus <7j 2 os 



a 2 there. In addition, of os 



where when observing the galaxy face-on (9 = 0), or in 
the case k = (no rotation). Finally, since the observed 
velocity dispersion is always measured within the aper- 
ture of the spectrograph, we integrate a 
apertures of radius R a : 



los over circular 



A/ P (i? a )a 2 (i? a ) ee / Z(x/,y/)af OB (x/,y/)dx>dy/. (11) 
JR'<R & 

In the following we will identify the quantity <7o appearing 
in the FP with cr a . 



2.3. The numerical code 

The solution of the Jeans equations and the projection of 
the density and of the various kinematical quantities are 
calculated with a grid-based numerical code developed by 
one of us (A.R.). The gravitational potential is obtained 
by solving the Poisson equation V 2 = 4nGp in spherical 
coordinates (r, ■&, <p), specializing the Londrillo & Messina 
(1990) spectral method to axisymmetric systems. The $ 
(co-latitude) dependence of p and <fi is described with stan- 
dard Legcndre polynomials, while the spherical radius is 
mapped as 



r ee tan 2 (m/2) , (0 < u < vr). 



(12) 



The resulting radial u-Laplace operator, when acting on 
Chebyshev polynomials T n (cos it) , originates a finite linear 
combination of such orthogonal functions, used to expand 
the density and the gravitational potential. Equations l@J 
and (0 are then integrated with a cubic spline algorithm 
on a logarithmic (R, z) grid, in order to sample the model 
central regions at high resolution (as required by the den- 
sity cusp and by the possible presence of a SMBH) . Since 
the projection coordinate system is tilted with respect to 
the natural coordinate system by an angle 9, the grid 
points of the two systems in general do not coincide, and 
so an interpolation and the los integrations are performed 
with a cubic spline algorithm. For obvious symmetry rea- 
sons, the computations are restricted to the first quadrant. 

The double precision C code is organized in a library of 
functions, and all the relevant intrinsic and projected fields 
of each model are computed in ~ 100 s on a 1.33GHz pro- 
cessor, for a 512 3 grid. The code accuracy has been tested 
by solving the Poisson, Jeans, and projected equations 
for the Miyamoto-Nagai (1975) models (see also Ciotti 
& Pellegrini 1996), for the spherical 7-models, and for the 
Ferrers ellipsoids: for these models all the relevant dynam- 
ical properties and, to some extent, even the projected 
fields are known analytically. In addition, we also com- 
pared our numerical results with the analytical (asymp- 
totic) models presented in the Appendix and in CB, ob- 
taining relative errors ft, 10~ 3 for the intrinsic dynamical 
properties and <10 -2 for the projected fields (also in the 
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central regions and in presence of a SMBH). Finally, for 
all models presented in the following section we verify nu- 
merically the projected virial theorem 

2n i n j K ij = -njW^ - n\W% 2 - n§W| 3 , (13) 

with Kij and Wij being the kinetic and potential energy 
tensors (e.g. Ciotti 2000, see also Eqs. [31], [A. 8], [A. 11], 
[A.12] in LC). 

3. General results 

Before studying in detail the effect of ordered rotation on 
the measured velocity dispersion in elliptical galaxies and 
its consequences on the FP, we present a few representa- 
tive models in order to illustrate the general properties of 
their projected kinematical fields (see also DG). In par- 
ticular, since ordered rotation enters in the definition of 
a 2 through the quantities v p and V p (Eq. ^U] and jlljl. 
and since these quantities are maximized for isotropic ro- 
tators seen edge-on, we restrict to this configuration in the 
following discussion. 

In Fig. n (top left panel) we show v p along the ma- 
jor axis for three different values of the galaxy density 
inner slope (7 = 0,1,1.5) and two different flattenings 
(q = 0.3, 0.8) in one-component isotropic models with- 
out central SMBH. Consistently with Eqs. I|A.5|) and ©, 
and with CB (Eq. [C.8]), in central regions v p vanishes in 
all the considered cases. The projected streaming velocity 
also vanishes independently of 7, as R — > 00, where the 
density profile oc to~ 4 . The major axis projected stream- 
ing velocity thus presents a rapid increase up to a max- 
imum, placed very near to the model circularized effec- 
tive radius, followed by a mild decrease for increasing ra- 
dial coordinate (note the logarithmic scale of the x-axis). 
Apparently, at fixed galaxy mass and scale-length, the 
main parameters determining the projected velocity v p 
are the galaxy flattening q, and the slope 7 of the central 
cusp. In Fig. n (top right panel) we show the aperture ve- 
locity dispersion cr a for the same models (with q = 0.3) in 
the left panel, either supported by tangential anisotropy 
(k — 0, solid line), or by ordered rotation (k = 1, dashed 
line). In general, cr a is larger for k = than for k = 1 
at any fixed aperture, and this is so because tangential 
anisotropy is maximal when k = 0, while the streaming 
motions (that are maximal for k = 1) contribute to <7i os 
with the term V p — v p which is small (see Eq. ^01)- This 
dependence of cr a on k is explicit in Eq. (29) in CB. In 
any case, also <7 a presents a maximum, which is near R c /8 
for the isotropic rotators, and near R e for non-rotating 
models. 

An important quantitative result of our models is that 
er a differs less than 15% between the non-rotating and 
maximally rotating cases, when R a <R e /8 since v p and 
V p vanish near the galaxy center. When integrating <7i os 
over a larger area (i? a ~ i? G ), however, the differences 
between <r a of the anisotropy-supported and the rotation- 
supported models increase for increasing 7, attaining a 
value of about 30% in the most extreme case (7 = 1.5). 



In Fig. n (bottom panels) we show the same quantities 
of the corresponding upper panels, for identical galaxy 
models at which we added a spherically symmetric (q^ = 
1) dark matter halo with Mh — 5Af» and Rh = 2R C . 
For the halo central cusp we adopted 7 = 1, in order to 
mimic the central behavior of the NFW density profile. 
Qualitatively the results are similar to those of the one- 
component models, the main effect being the substantial 
increase of <r a for large apertures in the k = case. Note 
also that in the central regions the effect of the adopted 
dark matter halo is negligible, as revealed by Eqs. (|A.4(1 - 
(IA.5|I with A/bh = 0. Analogous results are obtained also 
with halos with different mass or scale radius. 

Compared to the observations (e.g., Mehlert et al. 
2000), while the model rotational velocity profiles are 
quite realistic, the aperture velocity dispersion presents 
a central dip that seems at odd with what empirically in- 
ferred from the data (e.g., J0rgensen et al. 1995; Mehlert 
et al. 2003; Cappellari et al. 2005). This is an unrealistic, 
but very well known, feature of several dynamical mod- 
els with cuspy density distributions, which produce cen- 
trally vanishing velocity dispersion profiles (see Eq.[A.4]; 
e.g., Bailey & MacDonald 1981; Dehnen 1993; Tremaine 
et al. 1994; and Bertin et al. 2002 for a general discus- 
sion in the case of spherical symmetry). Such a central 
dip in the models often occurs at small radii, so that it 
is undetectable in the data or it nicely agrees with few 
observed cases (e.g., Graham et al. 1998; Emsellem et al. 
2004; Cappellari et al. 2005). This is also what we find 
for our isotropic rotators, but cannot apply to the case of 
non-rotating models, which show a peak in the los velocity 
dispersion profile well far from the galaxy center, near the 
effective radius. However, we stress again that here we are 
considering the most extreme (thus also possibly unrealis- 
tic) cases, in order to maximize the effects of ordered rota- 
tion, while these unobserved features become much weaker 
for milder flattenings, intermediate los inclinations, and 
non-zero rotational support. We therefore consider our re- 
sults reliable and conclude that the contribution of the 
rotational velocity to the observed central velocity disper- 
sion (i? a <i? e /8) is usually negligible for elliptical galaxies 
in the local Universe, also in presence of a cuspy dark mat- 
ter halo. This implies, that a systematic increase of rota- 
tional support with decreasing galaxy luminosity is not at 
the origin of the tilt of the FP of elliptical galaxies (see also 
Busarello et al. 1997, LC). However, our results also sug- 
gest that particular care should be used when constructing 
the FP of galaxies at high redshift (where spectroscopic 
apertures usually enclose a large fraction of the galaxy). In 
fact, the increasing difference of cr a between rotating and 
non rotating galaxies, with increasing i? a (dashed line in 
Fig. ^) would lead, if not properly taken into account, to 
underestimate the galaxy mass (when interpreting cr a by 
means of virial estimators based on spherical, not rotating 
models) up to 70% for low-mass and rotating objects (see 
also Bender, Burstein & Faber 1992, van Albada, Bertin & 
Stiavelli 1995). These aperture effects might so contribute 
(in part) to the observed decrease of the mass-to-light ra- 
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Fig. 1. Edgc-on, major axis projected streaming velocity v p (left panels), and aperture velocity dispersion a a (right 
panels) for one (upper panels) and two-component (lower panels) galaxy models. Projected streaming velocity is shown 
for isotropic models only, while er a for isotropic and fully tangcntially anisotropic (i.e. not rotating) models. Aperture 
radii corresponding to R c /8 and R c are marked as empty dots (o) and solid dots (•), respectively. In two-component 
models the dark matter halo is described by a spherical Hernquist (1990) model (7 = 1, q = 1 in Eq. [3]), with 
M h = 5 M* and R h = 2 R c . 



tio of low-mass galaxies at high redshift (e.g., di Serego 
Alighieri et al. 2005; Treu et al. 2005; van der Wei et al. 
2005). 

3.1. Effects of a central super-massive black hole 

It is now commonly accepted that SMBHs reside in the 
center of stellar spheroids, and that their mass Mbh scales 
almost linearly with the stellar mass of the parent galaxy, 
with Mbh/M* ps 10~ 3 (Magorrian et al. 1998). An order- 



of-magnitude estimate of the radius tbh of the "sphere of 
influence" (i.e., the region within which the gravitational 
effects of the SMBH on the stellar orbits are significant) is 
usually obtained as GMbh/tbh = ~ GM*/(2 r vir ), 

i.e., rBn/Re ~ 6 x 10~ 3 (where cr* is the one-dimensional 
stellar velocity dispersion, and r v i r /R c ~ 2.97 for a i? 1 / 4 
galaxy; see Ciotti 1991). Thus, SMBH effects should be de- 
tectable on scales of few parsecs or tens of parsecs only, in 
any case well within the standard central aperture of R c /8. 
However, while in galaxy models without SMBH v v and 
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0.01 0.1 1 0.01 0.1 1 

R I R c R J R c 

Fig. 2. Left panel: Intrinsic (solid line) and projected (dashed line) streaming velocities on the equatorial plane (z = 0) 
of a models with 7 = 1, q = 0.3, k = 1, and a SMBH with Mbh = 10~ 3 M* normalized to the corresponding quantity 
of the same model without the SMBH. Right panel: Intrinsic (solid line), line-of-sight (dashed line), and aperture 
velocity dispersion ratios (dotted line) for the same model. The aperture velocity dispersion is calculated for a circular 
aperture R a = R. The empty (o) and solid (•) dots correspond to R c /8 and R Cl respectively. 



Up vanish at the center for 7 < 2, in presence of a SMBH 
oc Wp H oc 1/yR, and in principle even an unresolved 
but luminosity weighted, central kinematical spike could 
produce detectable effects on cr a . The results obtained 
with our models are illustrated in Fig. (left panel), where 
we show the spatial and projected streaming velocity along 
the major axis of an isotropic 7 = 1 model with q = 0.3 
and with a central SMBH having M B h/M* = 10~ 3 . It 
can be seen that while v v increases by nearly one order 
of magnitude in the innermost regions, v p (the only one 
directly accessible to observations) is much less affected 
by the presence of the SMBH. Note that from Eq. ljA.51) . 
and CB (Eq. [C.8]), v® H /v v oc «™/«p oc R^ 3 ^ 2 . In 
all explored models we found that the SMBH influence is 
negligible for R>R c /8. 

In Fig. [21 (right panel) we show for the same model 
the intrinsic one-dimensional and los velocity disper- 
sions along the major axis, and the aperture velocity dis- 
persion as a function of i? a . Note that from Eq. (|A.4J| 
and Eqs. (C.3), (29)-(30) of CB we expect a^/an oc 
er^/erios oc erif H /(7 a oc R( J ~ 3 ^ 2 . In fact, our numerical 
results confirm this radial trend, and also show that the 
SMBH influence decreases along the sequence o~ z , <j\ os and 
er a at fixed radial coordinate. A comparison of the two 
panels in Fig. [21 shows that velocity dispersions are less 
affected than rotational velocities by the presence of a 
central SMBH. We thus conclude, that, in general, the 
presence of a SMBH can be neglected when using the 
measured "central" (i? a ~ Re/ 8) velocity dispersion, for 
the construction of the Faber-Jackson (1976), FP, and 
Mbh-<7 relations (Ferrarese & Merritt 2000, Gebhardt et 



al. 2000), as well as the u/u-ellipticity plane (Illingworth 
1977, Binney 1978). Accordingly, in the following Sections 
we will not consider a central black hole in the models. 

4. Contribution of rotational velocity to the FP 
thickness 

In a follow-up of LC we now investigate the contribution of 
projection and ordered rotation on the observed FP thick- 
ness. As anticipated in the Introduction, the only differ- 
ences with respect to the results of LC are possibly due to 
the new value of do , since the variations of the circularized 
effective radius R c (and of J e ) are independent of the spe- 
cific stellar density profile adopted (provided it is stratified 
on similar ellipsoids) . We start by illustrating how models 
"move" in the edge-on view of the FP as a function of 
their intrinsic (7, q, k) and observational (9) parameters. 
In Fig. |3 (left panel) we show the edge-on view of the FP 
of Coma cluster ellipticals as given by Eq. (solid line), 
and its rms scatter ~ 0.057 (dashed line, as obtained by 
JFK for galaxies with <tq > 100 km/s, and after correc- 
tion for measurement errors). The three families of mod- 
els superimposed are determined (from top to bottom) 
by (q, k) = (0.3, 0), (0.3, 1), and (0.6, 1). For simplicity, 
7 = 1 in all models, while three different aperture radii 
i? a are considered to measure cto = cr a , namely i? a = R c /8 
(dots), R c (squares), 10i? c (triangles). The cases corre- 
sponding to 8 = and i? a = R e are arbitrarily placed on 
the FP best-fit line, and the displacements correspond to 
6 increasing from to tt/2. According to Eq. (jHJ, when 
the los inclination changes from to 7r/2, R c decreases 
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(and / c increases), thus producing a vertical down-shifts 
towards the left of the representative models; as expected, 
displacements are smaller for rounder systems. The ad- 
ditional effect of rotational velocity produces horizontal 
displacements which depend mainly on the adopted R a . 
For k — 1 (isotropic rotators) , the variation of <7 a with the 
viewing angle 9 is only slightly dependent on the consid- 
ered spectroscopic aperture, so that the model displace- 
ments are almost parallel to each other, but not to the 
FP. On the contrary, the displacements of non-rotating 
anisotropic galaxies change considerably both in magni- 
tude and direction as function of R a . We have verified 
that these results are almost independent of the model 
central cusp 7. Thus, due to projection effects, galaxies 
move in directions which are not exactly parallel to the 
edge-on FP, but the entity of these displacements is small 
enough to always maintain the models within the observed 
FP thickness. A comparison with Fig. 5 in LC reveals re- 
markably similar behaviors. 

We now restrict the analysis to rotational effects as 
a function of density slope and galaxy flattening, and in 
Fig. |3 (right panel) we plot a selection of representative 
models. As in the left panel, we consider the edge-on view 
of the FP, and we describe how models "move" when vary- 
ing the relative amount of rotational and anisotropy sup- 
port. In particular we show the behavior of three families 
of models seen edge-on, with different central cusp slopes 
(7 = 0, 1, 1.5), and with different flattening q. Isotropic 
models (k = 1) are arbitrarily placed on the FP best-fit 
line. A reduction of the amount of rotational support to 
zero "moves" models toward the right, placing them at the 
intermediate points, corresponding to k = and q = 0.6. 
More flattened models move even further, until the right- 
most points, that represent models with k = and q = 0.3. 
These displacements are due to the increasing values of 
the central velocity dispersion ao as the rotational sup- 
port becomes negligible, and their amount increases with 
both the aperture radius and the galaxy flattening, while 
it is almost independent of the shape of the central cusp 7, 
in agreement with Fig. ^ Note that a preliminary analysis 
of this problem was given by CB, who adopted homeoidal 
expansion to investigate the behavior of the projected ve- 
locity dispersion in cuspy oblate models with a central 
SMBH. A comparison of our results with Fig. 2 of CB 
and results discussed therein, reveals very similar behav- 
iors. Note that the models showed in Fig. (right panel) 
are all observed edge-on, in order to maximize the rotation 
effect; when all the possible line-of-sight and the distribu- 
tion of intrinsic ellipticities are taken into account (as in 
the Monte-Carlo approach in LC) the statistical displace- 
ments would be essentially negligible. 

From the study of the models presented in this Section, 
and from their comparison with models in LC and CB we 
then confirm and extend the conclusions of LC to cuspy 
galaxies, i.e. that projection effects only marginally con- 
tribute to the FP thickness, 90% of it being due to vari- 
ations, from galaxy to galaxy, of their intrinsic physical 
properties. In particular, we proved that variations of ro- 



tational support from galaxy to galaxy only marginally 
contribute to the FP thickness. 

5. The v /cj-ellipticity plane 

With the aid of the developed models we finally study 
how the u/er-ellipticity plane (an important tool used to 
investigate to what extent the galaxy flattening is due to 
velocity dispersion anisotropy, e.g. see Illingworth 1977, 
Binney 1978) is populated by galaxies as a function of 
the adopted aperture i? a and of their streaming veloc- 
ity support. Traditionally, observational data are the peak 
projected streaming velocity t>p imax a "central" aperture 
velocity dispersion oq, and galaxy ellipticity. Data points 
are then compared with the family of curves obtained from 
the tensor virial theorem, where v 2 and a 2 are the mass 
weighted square streaming velocity and velocity disper- 
sion computed over the whole galaxy. As well known, the 
tensor virial theorem predictions are independent of the 
specific galaxy density profile as far as it is homeoidally 
stratified (Roberts 1962). For example, the thick solid line 

I in Fig. 0] is obtained from Eq. (4.95) of BT (see also 
Eq. [9] of Binney 2005), and corresponds to isotropic ro- 
tators seen edge-on. Galaxies whose representative points 
lie significantly below this curve are considered anisotropy 
supported. Unfortunately it is not straightforward to con- 
nect the virial quantities to the observable ones. For this 
reason, the tensor virial theorem curves are often "cor- 
rected" to take into account projection effects, and in 
general these corrections makes curve I flatter (e.g., see 
Eq. [4.5] and Figs. [8]-[9] in Evans & de Zeeuw 1994, and 
Figs. [4.5]-[4.6] in BT). The situation became even more 
complicate when Evans & de Zeeuw (1994) showed that 
v/a for their nearly isotropic "power-law" galaxy mod- 
els, not only was systematically below curve I of classical 
spheroids (and significantly so at large ellipticities), but 
also below the "corrected" locus. Of course, this behav- 
ior is not unexpected because, as discussed by Evans & 
de Zeeuw (1994), the density of their model is not strat- 
ified on ellipsoidal surfaces. In addition, their projected 
streaming velocity does not have a maximum, and its fidu- 
cial value must therefore be taken at some arbitrary dis- 
tance from the center. Interestingly, a very similar behav- 
ior is also shown by the isotropic models described in CB 
(Fig. [2]): in this case, due to the models scale- free na- 
ture, v p and cr a adopted to construct v/a were taken at 
the same radius. 

Fortunately, it is now possible to observationally map 
the rotation speed and the velocity dispersion over a sub- 
stantial fraction of a galaxy's image (e.g. the SAURON 
project, see de Zeeuw et al. 2002, Cappellari et al. 2005). 
Thus, one can use the sky-averaged quantities to define 
v/a, and predict its trend from the projected virial theo- 
rem (e.g. Binney 2005, see also Ciotti 1994, 2000). Curve 

II in Fig. 0] shows the corresponding locus of the edge-on 
isotropic rotators, as given by Binney (2005). 

Here we use our (isotropic, edge-on, and homeoidally 
stratified) models to study the effects of light profile cen- 
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Fig. 3. Left panel: Projection effects on the galaxy models in the space where the FP for Coma cluster galaxies is seen 
edge-on (solid line) with 1-rms scatter (dashed line) from JFK. The intrinsic model parameters are given by labels, 
and 7 = 1 for all models. For each model three spectroscopic apertures for er a = a are considered (Re/8: dots, R e : 
square, 10i? e : triangles). The los inclination increases from 9 = (face-on) to 9 = tt/2 (edge-on) along the arrows, and 
the cases 9 = 0, i? a = R c are arbitrarily placed on the FP best-fit line. Right panel: Maximum rotational effects as a 
function of intrinsic flattening and spectroscopic aperture i? a for three edge-on models (7 = 0, 1, 1.5). The isotropic 
rotators with q = 0.6 are arbitrarily placed on the FP while, on each horizontal line, the middle points correspond to 
models with k = and q = 0.6, while the rightmost points to models with k = and q = 0.3. 



tral cusps and different definitions of the projected stream- 
ing velocity, in determining the model position in the v/cr- 
e plane, while for simplicity we fix a = o~ a (R c /8). In par- 
ticular, we consider the cases v = v Pimax (the maximum 
projected velocity) and v = v p (R c /8), as representative 
of the measured rotation. The resulting trends are shown 
in Fig. 21 and it is clear how all the resulting curves lie 
below the (uncorrected) locus of isotropic, classical ellip- 
soids obtained from the tensor virial theorem (solid line I). 
Deviations with respect to the "classical" expectation are 
larger when v = v p (R c /8), a case analogous to those con- 
sidered by CB. Remarkably, our analysis shows that curve 
II (derived from the sky-averaged virial theorem) provides 
a much better estimate ofv/a (and so, of anisotropy) than 
the original one based on the (unprojected) tensor virial 
theorem, when a = <7 a (i? /8) and v = v P;max . Note also 
how the central cusp slope can be important, especially 
when using small aperture to determine v. 

Thus we conclude that particular care should be placed 
on the choice of the theoretical v/a locus to be used for 
comparison with the data, as a function of the quality and 
extension of the data themselves (see also DG). 

6. Discussion and conclusions 

In this paper we have quantified the contribution of or- 
dered rotation to the FP tilt and thickness. We adopted 
a class of oblate galaxy models with central cusp and ad- 



justable flattening and amount of rotational support. We 
also considered models with a dark matter halo and with 
a central super-massive black hole. The associated Jeans 
equations have been (numerically) solved under the as- 
sumption of an underlying two-integrals distribution func- 
tion. The main results can be summarized as follow: 

1. For the adopted family of models the effect of a central 
SMBH is stronger on streaming velocities than on ve- 
locity dispersions. However, for mass ratios of the order 
of those predicted by the Magorrian relation, the effect 
is in general negligible, and the presence of a central 
SMBH can be neglected when studying relations such 
as the FP. 

2. In agreement with previous results of LC and CB, also 
for the adopted class of cuspy models the contribution 
of streaming motions to the observed velocity disper- 
sion is negligible when small/medium apertures (< R e ) 
are used for the spectroscopic observations. This im- 
plies that a systematic decrease of rotational support 
with increasing luminosity is not at the origin of the 
tilt of the FP of elliptical galaxies at low-redshift. 

3. In general, the aperture velocity dispersion measured 
at R a >R c is larger for not rotating models than 
for rotationally supported models (and the difference 
increases significantly when a dark matter halo is 
present). This must be taken properly into account 
when studying the FP, or estimating the mass-to-light 
ratios at high redshift. 
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Fig. 4. Rotational parameter for isotropic, edge-on models 
as a function of their ellipticity. The solid thick curve I is 
the locus of isotropic, classical ellipsoids as obtained from 
the tensor virial theorem (Binncy 1978; see also Eq.[4.95] 
in BT) , while solid line II is derived from the sky-averaged 
tensor virial theorem (Eq. [26] of Binney 2005, with 5 = 
and a = 0) . The upper set of curves represents our models 
when v p = v Pima x, while the lower set of curves represents 
our models when v = v p (R c /8). In all cases we adopt 

a = C7 a (R c /8). 



4. When observed from different los inclinations, mod- 
els move in the (logi? e , log/ e , logco) space along di- 
rections that are not parallel to the edge-on view of 
the FP, thus confirming that projection effects do 
contribute to the observed FP scatter. The models 
displacement depends mainly on their intrinsic flat- 
tening (being larger for more flattened systems), on 
the amount of rotational support, and on the aper- 
ture used to measure <jq. These effects however, when 
weighted over the distribution of cllipticities and los 
inclinations, would produce statistical displacements 
well within the FP thickness, in agreement with the 
conclusions of LC. 

5. Finally, as a by-product of the present investigation 
we studied the position of the models in the vja- 
cllipticity plane, a tool used to discriminate between 
rotationally and velocity dispersion supported galax- 
ies. We found that, while our models remain well below 
the (uncorrected) locus of classical ellipsoids as deter- 
mined from the tensor virial theorem, their position 
is well predicted by sky-averaged projected virial the- 
orem (Binney 2005) if the peak rotational velocity is 
used. 
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Appendix A: A simple, analytical, cuspy and 
axisymmetric galaxy model with dark matter 
halo and central SMBH 

We present here the internal kinematical fields of a family 
of power-law galaxy models (constructed with a homeoidal 
expansion following CB) with dark matter halo and cen- 
tral black hole. The seed distributions for stars and dark 
matter are p„ = p*o/ m 2 and ph = pho/m^, where 



m 



2 _ 



= R 2 + z 2 /(l - e) 



mi 



Ph 

it 2 



z : /(l - rj) and, 



without loss of generality, the scale-length of both distri- 
butions is R Cl with R = R/R c and z = z/R c . According 
to CB (Eq. [A. 2]), the (unconstrained) density and po- 
tential expansion are performed for < 1, dr] < 1 (with 
7, 5 > 1). Thus, from CB (Eq. [26]), the model stellar 
density distribution and the total potential are 
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P* = 



je jeR 



py+2 
P 

- ) 

r 



(A.l) 
(A.2) 



where the density is normalized to p *, the gravitational 
potentials to AirGp^R 2 , p = Mq^/Attp^oRc, an d finally 
Tl = Pim/p*o- The explicit expression for and ah h is 
given in CB (Eq. [27]), and the model (dimensionless) cir- 
cular velocity is 
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The Jeans equations for the stellar component are inte- 
grated following CB[Eqs.(14)-(15)], retaining only first or- 
der terms in the flattenings e and m 
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